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Abstract. We discuss the the notion of a partial dynamical symmetry (PDS), for which a 
prescribed symmetry is obeyed by only a subset of solvable eigenstates, while other eigenstates 
are strongly mixed. We present an explicit construction of Hamiltonians with this property, 
including higher-order terms, and portray their significance for spectroscopy and shape-phase 
transitions in nuclei. The occurrence of both a single PDS, relevant to stable structures, and of 
several PDSs, relevant to coexistence phenomena, are considered. 



1. Introduction 

Models based on spectrum generating algebras form a convenient framework to examine un- 
derlying symmetries in dynamical systems, and have been used extensively in diverse areas of 
physics pQ. Notable examples in nuclear physics are Wigner's spin-isospin SU(4) supermulti- 
plets [2], SU(2) single- j pairing [3], Elliott's SU(3) model [I], symplectic model [5], Ginocchio's 
monopole and quadrupole pairing models [6], interacting boson models (IBM) for even-even 
nuclei [7] and boson- fermion models (IBFM) for odd-mass nuclei [S|. Similar algebraic tech- 
niques have proven to be useful in the structure of molecules [SJ[TU] and of hadrons [llj . In 
such models the Hamiltonian is expanded in elements of a Lie algebra, (Go), called the spec- 
trum generating algebra. A dynamical symmetry occurs if the Hamiltonian can be written in 
terms of the Casimir operators of a chain of nested algebras, Go D Gi D . . . D G n [12J. The 
following properties are then observed, (i) All states are solvable and analytic expressions are 
available for energies and other observables. (ii) All states are classified by quantum numbers, 
|ao,ai, . . . ,a n ), which are the labels of the irreducible representations (irreps) of the algebras 
in the chain, (iii) The structure of wave functions is completely dictated by symmetry and is 
independent of the Hamiltonian's parameters. 

A dynamical symmetry provides clarifying insights into complex dynamics and its merits 
are self-evident. However, in most applications to realistic systems, the predictions of an exact 
dynamical symmetry are rarely fulfilled and one is compelled to break it. The breaking of the 
symmetry is required for a number of reasons. First, one often finds that the assumed symmetry 
is not obeyed uniformly, i.e., is fulfilled by only some of the states but not by others. Certain 
degeneracies implied by the assumed symmetry are not always realized, (e.g., axially deformed 
nuclei rarely fulfill the IBM SU(3) requirement of degenerate /3 and 7 bands [7]). Secondly, 
forcing the Hamiltonian to be invariant under a symmetry group may impose constraints which 
are too severe and incompatible with well-known features of the dynamics (e.g., the models 
of [6] require degenerate single- nucleon energies). Thirdly, in describing systems in-between two 
different structural phases, e.g., spherical and deformed nuclei, the Hamiltonian by necessity 
mixes terms with different symmetry character. In the models mentioned above, the required 



Table 1. Generators, Casimir operators, Ck(G), of order k = 1,2,3 and their eigenvalues for algebras G 
in the IBM. Here n s = sh, h d = s/5U(°\ N = h s + h d , L m = y/WuS?, Q m = - ^U%\ 
11$ = dls + S trf m , fi^ = i{d\ n s - S trf m ), U%> = (dt d)#, where d m = (-l) m d. m . 
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symmetry breaking is achieved by including in the Hamiltonian terms associated with (two or 
more) different sub-algebra chains of the parent spectrum generating algebra. In general, under 
such circumstances, solvability is lost, there are no remaining non-trivial conserved quantum 
numbers and all eigenstates are expected to be mixed. A partial dynamical symmetry (PDS) [13] 
corresponds to a particular symmetry breaking for which some (but not all) of the virtues of 
a dynamical symmetry are retained. The essential idea is to relax the stringent conditions of 
complete solvability so that the properties (i)-(iii) are only partially satisfied. It is then possible 
to identify several types of partial dynamical symmetries. PDS of type I corresponds to a 
situation where some of the states have all the dynamical symmetry. In this case, properties 
(i)-(iii) are fulfilled exactly, but by only a subset of states. PDS of type II corresponds to a 
situation for which all the states preserve part of the dynamical symmetry. In this case, there 
are no analytic solutions, yet selected quantum numbers (of the conserved symmetries) are 
retained. PDS of type III has a hybrid character, for which some of the states preserve part of 
the dynamical symmetry. 

In what follows we discuss algorithms for constructing Hamiltonians with partial dynamical 
symmetries and demonstrate their relevance to quantum systems. For that purpose, we employ 
the interacting boson model (IBM) [7], widely used in the description of low-lying quadrupole 
collective states in nuclei in terms of interacting monopole (s) and quadrupole (d) bosons 
representing valence nucleon pairs. The bilinear combinations {s^s, s^d m , dlnS, d\nd m '} span 
a U(6) algebra, which serves as the spectrum generating algebra. The IBM Hamiltonian is 
expanded in terms of these generators and consists of Hermitian, rotational-scalar interactions 
which conserve the total number of s- and d- bosons, N = n s + hd = s^s + Yl m d\nd m - Three 
dynamical symmetry limits occur in the model with leading subalgebras U(5), SU(3), and 
0(6), corresponding to typical collective spectra observed in nuclei, vibrational, rotational, 
and 7-unstable, respectively. Relevant information on these algebras is collected in Table [TJ 



A geometric visualization of the model is obtained by an energy surface 



E N (p^) = <j3, r ,N\A\p,r,N) , (l) 

defined by the expectation value of the Hamiltonian in the coherent (intrinsic) state |14yi5| 

|/3, 7; AO = (N\)- l / 2 (bl) N \0) , (2a) 
b\ = (l + /3 2 )~ 1 / 2 [/3cos74 + ^sin 7 (4 + dl 2 )/^ / 2 + st ] • ( 2b ) 



Here (yS, 7) are quadrupole shape parameters whose values, (/3o>7o)> at the global minimum 
of £Jjv(/3,7) define the equilibrium shape for a given Hamiltonian. The shape can be spherical 
[fi = 0) or deformed (/3 > 0) with 7 = (prolate), 7 = tt/3 (oblate), < 7 < 7r/3 (triaxial), or 
7-independent. The equilibrium deformations associated with the dynamical symmetry limits 
are fa = for U(5), (/3 = \/2,7o = 0) for SU(3) and (/3 = 1, 70 arbitrary) for 0(6). 

2. Construction of Hamiltonians with partial dynamical symmetries 

PDS of type I corresponds to a situation for which the defining properties of a dynamical 
symmetry (DS), namely, solvability, good quantum numbers, and symmetry-dictated structure 
are fulfilled exactly, but by only a subset of states. An algorithm for constructing Hamiltonians 
with PDS has been developed in |16j and further elaborated in [T7]. The analysis starts from 
the chain of nested algebras 

Gdyn ^ G D D G S y m 

III (3) 
[h] <£> A 

where, below each algebra, its associated labels of irreps are given. Eq. ([3]) implies that Gd yn is 
the dynamical (spectrum generating) algebra of the system such that operators of all physical 
observables can be written in terms of its generators; a single irrep of Gd yn contains all states 
of relevance in the problem. In contrast, G sym is the symmetry algebra and a single of its irreps 
contains states that are degenerate in energy. Assuming, for simplicity, that particle number 
is conserved, then all states, and hence the representation [h], can then be assigned a definite 
particle number N. For N identical particles the representation [h] of the dynamical algebra 
Gdyn is either symmetric [N] (bosons) or antisymmetric [1^] (fermions) and will be denoted, in 
both cases, as [hj^]. The occurrence of a DS of the type signifies that the Hamiltonian is 
written in terms of the Casimir operators of the algebras in the chain, 

H D s = Y, a G C ( G ) , (4) 

G 

and its eigenstates can be labeled as \[h]\r](T,) . . . A); additional labels (indicated by ...) are 
suppressed in the following. The eigenvalues of the Casimir operators in these basis states 
determine the eigenenergies £ , £>5([/itv](S)A) of Hds- Likewise, operators can be classified 
according to their tensor character under Q as T^^^x- 

Of specific interest in the construction of a PDS associated with the reduction (|3|), are the 
n-particle annihilation operators T which satisfy the property 

r [M(CT)A |[M(s )A) = o, (5) 

for all possible values of A contained in a given irrep (So) of G. Equivalently, this condition can 
be phrased in terms of the action on a lowest weight (LW) state of the G-irrep (So), 



f [hn]{a)x \LW; [M(£o» = 0, 



(6) 



from which states of good A can be obtained by projection. Any re-body, number-conserving 
normal-ordered interaction written in terms of these annihilation operators and their Hermitian 
conjugates (which transform as the corresponding conjugate irreps), 

H> = Y,A a pfi% , (7) 

a,/3 

has a partial G-symmetry. This comes about since for arbitrary coefficients, A a p, H' is not a 
G-scalar, hence most of its eigenstates will be a mixture of irreps of G, yet relation ([5]) ensures 
that a subset of its eigenstates |[/ijv](So)A), are solvable and have good quantum numbers under 
the chain ([3]). An Hamiltonian with partial dynamical symmetry is obtained by adding to H' 
the dynamical symmetry Hamiltonian, Hds ©> still preserving the solvability of states with 
<S) = (So), 

Hpds = Hds + H . (8) 

If the operators Iry/jU span the entire irrep (a) of G, then the annihilation condition ^ is 
satisfied for all A-states in (So), if none of the G irreps (S) contained in the Gd yn irrep [Ii^-n] 
belongs to the G Kronecker product (a) x (So). So the problem of finding interactions that 
preserve solvability for part of the states Q is reduced to carrying out a Kronecker product. 
The arguments for choosing the special irrep (S) = (So) in Eq. ([5]), which contains the solvable 
states, are based on physical grounds. A frequently encountered choice is the irrep which contains 
the ground state of the system. The above algorithm is applicable to any semisimple group. 

PDS of type II corresponds to a situation for which all the states of the system preserve 
part of the dynamical symmetry, Go D G\ D G2 D ■ ■ ■ D G n . In this case, there are no 
analytic solutions, yet selected quantum numbers (of the conserved symmetries) are retained. 
This occurs, for example, when the Hamiltonian contains interaction terms from two different 
chains with a common symmetry subalgebra [18J, e.g., 

G D I ^ I D G 2 D . . . D G n . (9) 

If G\ and G[ are incompatible, i.e., do not commute, then their irreps are mixed in the 
eigenstates of the Hamiltonian. On the other hand, since G2 and its subalgebras are common 
to both chains, then the labels of their irreps remain as good quantum numbers. 

An alternative situation where PDS of type II occurs is when the Hamiltonian preserves only 
some of the symmetries Gi in the DS chain and only their irreps are unmixed. A systematic 
procedure for identifying interactions with such property was proposed in [15] . Let G\ D G2 D 
G3 be a set of nested algebras which may occur anywhere in the chain, in-between the spectrum 
generating algebra Go and the invariant symmetry algebra G n . The procedure is based on 
writing the Hamiltonian in terms of generators, gi, of G\, which do not belong to its subalgebra 
G2. By construction, such Hamiltonian preserves the G\ symmetry but, in general, not the G2 
symmetry, and hence will have the G\ labels as good quantum numbers but will mix different 
irreps of G%. The Hamiltonians can still conserve the G3 labels e.g., by choosing it to be a scalar 
of G3. The procedure involves the identification of the tensor character under G% and G3 of the 
operators gi and their products, gigj . . .gt- The Hamiltonians obtained in this manner belong 
to the integrity basis of G3-scalar operators in the enveloping algebra of G\ and, hence, their 
existence is correlated with their order. 

PDS of type III combines properties of both PDS of type I and II. Such a generalized PDS [20J 
has a hybrid character, for which part of the states of the system under study preserve part of the 
dynamical symmetry. In relation to the dynamical symmetry chain of Eq. ([3]), with associated 



basis, | [/ijv](S)A), this can be accomplished by relaxing the condition of Eq. ([5]), so that it holds 
only for selected states A contained in a given irrep (So) of G and/or selected (combinations of) 
components A of the tensor TJ/, n ]U)x- Under such circumstances, let G' 7^ G sym be a subalgebra 
of G in the aforementioned chain, G D G' . In general, the Hamiltonians, constructed from these 
tensors, in the manner shown in Eq. ©, are not invariant under G nor G' . Nevertheless, they 
do posses the subset of solvable states, |[/iat](Eo)A), with good G-symmetry (Eo) (which now 
span only part of the corresponding G-irrep), while other states are mixed. At the same time, 
the symmetry associated with the subalgebra G', is broken in all states (including the solvable 
ones). Thus, part of the eigenstates preserve part of the symmetry. These are precisely the 
requirements of PDS of type III. 



3. SU(3) partial dynamical symmetry 

The SU(3) DS chain of the IBM and related quantum numbers are given by [7] 

U(6) D SU(3) D 0(3) 
I I I (10) 

[N] (A, /i) K L 

For a given U(6) irrep [N], the allowed SU(3) irreps are (A, y) = (2N — 4k — 6m, 2k) with k, m 
non-negative integers, such that, A, y > 0. The multiplicity label K is needed for complete 
classification and corresponds geometrically to the projection of the angular momentum on 
the symmetry axis. The values of L contained in a given SU(3) irrep (A, /i), are obtained 
from the known SU(3) D 0(3) reduction. The states | [N] (A, fi)K L) form the (non-orthogonal) 
Elliott basis [I] and the Vergados basis |[iV](A, y)xL) [7] is obtained from it by a standard 
orthogonalization procedure. The two bases coincide in the large-N limit and both are eigenstates 
of a Hamiltonian with SU(3) DS. The latter, for one- and two-body interactions, can be 
transcribed in the form 



Hds = h 2 



G 2 (SU(3)) + 2N(2N + 3) +GG 2 (0(3)), (11) 



where G 2 (G) is the quadratic Casimir operator of G, as defined in Table [TJ The spectrum of 
Hds is completely solvable with eigenenergies 

E DS = h 2 6 [2N(k + 2m) - k(2k - 1) - 3m(2m - 1) - 6km] + CL(L + 1) , (12) 

and (A, fi) = {2N — 4k — 6m, 2k). The spectrum resembles that of an axially-deformed 
rotovibrator and the corresponding eigenstates are arranged in SU(3) multiplets. In a given 
SU(3) irrep (A,/x), each K- value is associated with a rotational band and states with the same 
L, in different X-bands, are degenerate. The lowest SU(3) irrep is (2N, 0), which describes the 
ground band g[K = 0) of a prolate deformed nucleus. The first excited SU(3) irrep (2A^ — 4, 2) 
contains both the j3{K = 0) and 7 (if = 2) bands. States in these bands with the same angular 
momentum are degenerate. This /3-7 degeneracy is a characteristic feature of the SU(3) limit of 
the IBM which, however, is not commonly observed. In most deformed nuclei the (3 band lies 
above the 7 band. In the IBM framework, with at most two-body interactions, one is therefore 
compelled to break SU(3) in order to conform with the experimental data. 

The construction of Hamiltonians with SU(3)-PDS of type I is based on identification of 
n-boson operators which annihilate all states in a given SU(3) irrep (A, fi), chosen here to be 
the ground band irrep (2N, 0). For that purpose, we consider the following two-boson SU(3) 
tensors, ^Lin^ M )x-£m> w ith n = 2, (A, fi) = (0, 2) and angular momentum £ = 0, 2 

■^[2](0,2)0;00 ^ ^0 = ■ S - 2(s 1 ") 2 , (13a) 

4](0,2)0; 2m « PL = ^ + V7(dU^ . (13b) 



The corresponding Hermitian conjugate boson-pair annihilation operators, Po an d P 2m , 
transform as (2,0) under SU(3), and satisfy 



P \[N](2N,0)K = 0,L) = 0, 
P 2m \[N](2N,0)K =0,L) = 0, £ = 0,2,4,...,2iV . (14) 

Equivalently, these operators satisfy 

Po|/3 = V2, 7 = 0;iV) = 0, 
P 2m \P = V2, 1 = 0;N) = 0, (15) 

where |/3 = v2, 7 = 0; JV), is the condensate of Eq. ([2]) with the SU(3) equilibrium deformations. 
It is the lowest-weight state in the SU(3) irrep (A, fi) = (2N, 0) and serves as an intrinsic state 
for the SU(3) ground band. The rotational members of the band | [N] (2N, 0)K = 0, L) , Eq. f)14[) , 
are obtained from it by 0(3) projection, and span the entire SU(3) irrep (A,/i) = (2iV, 0). The 
relations in Eqs. (fl~4"|) - (fl~5"]) follow from the fact that the action of the operators Pi m leads to 
a state with N — 2 bosons in the U(6) irrep [N — 2], which does not contain the SU(3) irreps 
obtained from the product (2,0) x (27V, 0) = (2iV + 2,0) (2N, 1) © (2iV - 2,2). In addition, 
Pq satisfies 

P \[N](2N -4k,2k)K = 2k,L) =0 , L = K, K + 1, . . . , (2N - 2k) . (16) 

For k > the indicated L-states span only part of the SU(3) irreps (A, fx) = (2N — 4k, 2k) and 
form the rotational members of excited 7 (K = 2k) bands. This result follows from the fact that 
Pq annihilates the intrinsic states of these bands, \ j k (K = 2k)) oc (Pj 2 ) fc |/3 = \/2, 7 = 0;N — 2k). 

Following the general algorithm, a two-body Hamiltonian with partial SU(3) symmetry can 
now be constructed as in Eq. ([7|), H' = hoP^Po + /12-pj ' where p2 m = (— ) m P2- m - For 
h 2 = ho, this Hamiltonian is an SU(3) scalar, while for ho = —5h 2 , it transforms as a (2,2) 
SU(3) tensor component. The scalar part is related to the quadratic Casimir operator of SU(3) 

O2 = P^Po + P 2 ■ Pi = -C 2 (SV(3)) + 2N(2N + 3) , (17) 

and is simply the first term in Hds-, Eq. (fTTj) . In accord with Eq. ([8]), the two-body SU(3)-PDS 
Hamiltonian is thus given by 

Hpds = Hds + tjP^Pq. (18) 

The PqPq term is not diagonal in the SU(3) chain, however, Eqs. (fT4"ft - (fl~6"|) ensure that Hpds 
retains selected solvable states with good SU(3) symmetry. Specifically, the solvable states are 
members of the ground g(K = 0) 

\N,(2N,0)K = 0,L) £ = 0,2,4, . . . ,2N (19a) 
E PDS = CL(L + 1) (19b) 



and j k (K = 2k) bands 



\N,(2N -4k,2k)K = 2k,L) L = K,K + 1, K + 2, . . . , (2N - 2k) (20a) 
£ , PDS = /i2 6/c(A^-2)(2iV-2A; + l) + C£(£ + l) k > . (20b) 



The remaining eigenstates of Hpds do not preserve SU(3) and, therefore, get mixed. 
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Figure 1. Spectra of 168 Er (N — 16). Experimental energies (EXP) are compared with IBM 
calculations in an exact SU(3) dynamical symmetry [SU(3)], in a broken SU(3) symmetry (WCD) and 
in a partial dynamical SU(3) symmetry (PDS). The latter employs the Hamiltonian of Eq. (fT5|) . with 
h 2 = 4, r) = 4, C = 13 keV gT]. 



The empirical spectrum of 168 Er is shown in Fig. [T] and compared with SU(3)-DS, SU(3)- 
PDS and broken SU(3) calculations [2T]. The SU(3)-PDS spectrum shows an improvement over 
the schematic, exact SU(3) dynamical symmetry description, since the ^-7 degeneracy is lifted. 



Table 2. B(E2) branching ratios from states in the 7 band in 168 Er. The column EXP lists the 
experimental ratios, PDS is the SU(3) partial dynamical symmetry calculation and WCD is a broken 
SU(3) calculation gj. 
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The quality of the calculated PDS spectrum is similar to that obtained in the broken-SU(3) 
calculation, however, in the former the ground g(K = Oi) and j(K = 2i) bands remain solvable 
with good SU(3) symmetry, (A,/x) = (2N, 0) and (2N — 4,2) respectively. At the same time, the 
excited K = 0^ band involves about 13% SU(3) admixtures into the dominant (2N — 4, 2) irrep. 
Since the wave functions of the solvable states (|19p -(|20 p are known, one can obtain analytic 
expressions for matrix elements of observables between them. For example, the most general 
one-body E(2) operator reads T(E2) = aQ + pll^ 2 ), in the notation of Table [TJ Since Q is 
an SU(3) generator, it cannot connect different SU(3) irreps, hence only H^ 2 \ which is a (2,2) 
SU(3) tensor, contributes to 7 — > g transitions. Accordingly, the calculated B(E2) ratios for 
7—^5 transitions involve ratios of known SU(3) isoscalar factors and lead to parameter-free 
predictions. The latter, as shown in Table [51 are in excellent agreement with experiment, thus 
confirming the relevance of SU(3)-PDS to the spectroscopy of 168 Er [21] . 

The construction of SU(3)-PDS Hamiltonians with higher-order terms follows the general 
algorithm and is based on identification of n-boson operators which annihilate all states in the 
irrep (2N,0). For n = 3, we consider the following SU(3) tensors, B\ n \(\ ^-t^ 

^[3](2,2)0;00 a Wq = hP^S^ — P 2 ■ d) , ^[ 3 ]( 2j 2)2;2m ^ ^2m = ^O^ln + ^P^S^ , 

B\ m ^2 m ^ V L = ZPUl-Pl m J , Bl m2)2 . em <xW} m = (P^ £ = 3A (21a) 
£[3](0,0)0;00 « At=P t s t + pt. d t. ( 21b) 

The operators wj, w\ m , v} m , W^ m , w\ m in Eq. ([21a]) span the irrep (A,//) = (2,2), while A f 
of Eq. (|21bp transforms as (A,//) = (0,0). The latter SU(3)-scalar operator is related to the 
cubic and quadratic Casimir operators of SU(3), defined in Table [H 

2A 1 A = C 3 (SU(3)) - 2iV(4iV + 3)(2A> + 3) + 3(2A + 3)£ 2 , (22) 

where 9 2 is given in Eq. (j!7p . In the presence of two- and three-body terms, the dynamical- 
symmetry Hamiltonian and eigenenergies for states with (A, (J,) = (2N — 4k — 6m, 2k), read 

H DS = hi At A + h 2 6 2 + C L ■ L , (23a) 
E DS = +h 1 54m{N(2k + 2m + l)-k(2k-l)-(2m-l)(2m + l)-6km] 

+h 2 6 [ 2N(k + 2m) - k(2k - 1) - 3m(2m - 1) - 6km } + C L(L + 1) . (23b) 



The SU(3)-PDS Hamiltonian has the following SU(3)-tensor expansion 

Hpds = Hds + V PqPo + ai W^W + a 2 (wjA + At Wo) + a 3 W\ ■ W 2 

+a 4 V% -V 2 + a 5 (w\ ■ V 2 + V ] 2 ■ W 2 ^j + a 6 w\ ■ W 3 + a 7 w\ ■ W 4 . (24) 

The relations of Eq. ([T4"]) ensure that the operators Wg m , V 2m and A of Eq. ([2T]) annihilate the 
states of the SU(3) ground band g(K = 0). The solvable eigenstates and eigenenergies of Hpds 
are those shown in Eq. (j!9p . In addition, the operator A (I21bp annihilates all states in the irreps 
{2N -4k,2k) 



A\[N]{2N - Ak,2k)K;L) = 



(25) 



This property follows from the fact that the U(6) irrep [N — 3] does not contain SU(3) irreps 
obtained from the product (0, 0) x (2N — 4k, 2k). Using Eqs. f)16|) and ()25p . we can identify the 
following sub-class of SU(3)-PDS Hamiltonians 

Hpds-i = H D s + h 3 P£Po + h 4 P^sP + h 5 (At s p + 4 stA ) , ( 26 ) 

with additional solvable states which are the members of the g(K = 0) and ^ k (K = 2k) bands 
listed in Eqs. (fT9]) - (|20|) . 

A second sub-class of SU(3)-PDS Hamiltonian corresponds to the choice 

Hpds-2 = H DS + h 6 W\ -W 2 + h 7 Wl ■ W 3 . (27) 

Hpds-2 has a solvable ground band g(K = 0), Eq. (fl~9j) . and a solvable (3{K = 0) band 

\N,(2N -4,2)K = 0,L) L = 0,2,4, ... , (2N — 4) (28a) 
Epds-2 = h 2 6(N - 2)(2N - 1) + CL(L + 1) . (28b) 

This result follows from the fact that Wim an d W 3m annihilate the intrinsic state of this band, 
|/3 (if = 0)) oc (x/2-Po ~ P 2,o)\P = V2, 7 = 0; N - 2), as well as the projected Elliott basis states 

W em \[N](2N -4,2),K = 0,L) = £ = 2,3. (29) 

Three-body terms allow an additional solvable symmetry-conserving operator 

Cl = —4V3Q • (L x L) (2) . (30) 

This operator is constructed from SU(3) generators, hence is diagonal in (A, fx). It breaks, 
however, the aforementioned if-degeneracy. A well defined procedure exists for obtaining the 
eigenstates of (l and corresponding eigenvalues (Cl) [22,23j. For example, for the irreps (A, 0) 
and (A, 2) with A even, we have 

(A,0)ir = 0, L = 0,2,4,...,A: (Cl) = (2A + 3)L(L + 1) , (31a) 

(X,2)K = 2, L = 3,5,7,... ,A + l,A + 2 : {Cl) = (2A + 5)[L(L + 1) - 12] , (31b) 

(A,2)A" = 0, L = 0: (Cl) = , (31c) 

(A, 2) K = 0,2, L = 2,4,6, ...,A : 

(Cl) = (2A + 5) [L(L + 1) - 6] ± 6 V / (2A + 5) 2 + L(L + 1)(L - 1)(L + 2) . (31d) 

Several works have examined the influence of the symmetry-conserving operator fi (|30p on 
nuclear spectra, within the IBM [23-25] and the symplectic shell model |26p27j . It is interesting 
to note that the operator Cl can be expressed in terms of Hpds-i and Hpds-2 (123) as 

Cl = -2(2N - 1% + 2A f A + (4A> + 3)L • L 

+6(A> - l)P f P - 2 (kUP + P ] Q s ] A) + 2 W 2 f • W 2 + 4 W 3 f • #3 • (32) 

4. 0(6) partial dynamical symmetry 

The 0(6) DS chain of the IBM and related quantum numbers are given by [7] 

U(6) D 0(6) D 0(5) D 0(3) 
I I I 4 (33) 

[N] (S) (r) n A L 



For a given U(6) irrep [N], the allowed 0(6) and 0(5) irreps are £ = N, N — 2, ... or 1, and 
t = 0, 1, . . . £, respectively. The values of L contained in the 0(5) irrep (r) are obtained from 
theknownO(5) D 0(3) reduction and ba is a multiplicity label. The eigenstates |[iV](E)(r)nA-^) 
are obtained with a Hamiltonian with 0(6) DS which, for one- and two-body interactions, can 
be transcribed in the form 



Hds = ho 



C 2 (0(6)) + iV(iV + 4) + 5 C 2 (0(5)) + C C 2 (0(3)) . (34) 



Here the quadratic Casimir operators, C*2(G), are defined in Table [TJ The spectrum of Hbs is 
completely solvable with eigenenergies 

Pds = 4/i (iV-t; + 2)t) + BT(r + 3)+CL(I + l). (35) 

The spectrum resembles that of a 7-unstable deformed rotovibrator, where states are arranged 
in bands with 0(6) quantum number £ = N — 2v, (v = 0, 1, 2, . . .). The ground band (v = 0) 
corresponds to the 0(6) irrep with £ = N. The 0(5) and 0(3) terms in Hbs (|34p . govern 
the in-band rotational splitting. The lowest members in each band have quantum numbers 
(r = 0, L = 0), (r = 1, L = 2) and (r = 2, L = 2, 4). 

The construction of Hamiltonians with 0(6)-PDS of type I is based on identification of n- 
boson operators which annihilate all states in a given 0(6) irrep, (£), chosen here to be the 
ground band irrep (£) = (N). For that purpose, a relevant operator to consider is 

4]<0>(0)0;00 « P W • dt _ ( S t)2 . (36) 

The corresponding Hermitian conjugate boson-pair annihilation operator, Pq, transforms also 
as (E) = (0) under 0(6) and satisfies 

P \{N}(N)(T)n A L)=0, r = 0,l,2,...,JV. (37) 

Equivalently, this operator satisfies 

Pol/? = 1,7; AO = 0, (38) 

where |/3 = 1,7; N), is the condensate of Eq. (J2J with the 0(6) equilibrium deformations. It is 
the lowest- weight state in the 0(6) irrep (£) = (N) and serves as an intrinsic state for the 0(6) 
ground band. The rotational members of the band, \ [N]{N) (t)tl/\L) , Eq. (|37p . are obtained from 
it by 0(5) projection, and span the entire 0(6) irrep (E) = (N). The relations in Eqs. ([37|) - ([38|) 
follow from the fact that the action of the operator Pq leads to a state with N — 2 bosons in 
the U(6) irrep [N — 2], which does not contain the 0(6) irrep (N) obtained from the product of 
(O)x(iV). 

Since both Pq and Po (|36p are 0(6) scalars, they give rise to the following interaction 

PjPo = -Cb( 6 ) + N(N + 4) , (39) 

which is simply the 0(6) term in Pds> Eq. (|34p . with an exact 0(6) symmetry. Thus, in this case, 
unlike the situation encountered with SU(3)-PDS, the algorithm does not yield an 0(6)-PDS 
of type I with two-body interactions. In the IBM framework, an Hamiltonian with a genuine 
0(6)-PDS of this class, requires higher-order terms. 

Focusing on three-body interactions with 0(6)-PDS, one considers the following two three- 
boson 0(6) tensors, Bt-,,i,*^ m , with n = 3, a = 1 and t = 0, 2, 

"^[3](1)(0)0;00 X P S ^ ' ( 40a ) 
^[3](l)(l)0;2m X Po4n • (40b) 
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Figure 2. Observed spectrum of 196 Pt compared with the calculated spectra of Hus (|34[) . with 
0(6) dynamical symmetry (DS), and of Hpus JUT) with 0(6) partial dynamical symmetry (PDS). The 
parameters in H US (#pds) are h = 43.6 (30.7), B = 44.0 (44.0), C = 17.9 (17.9), and rj = (8.7) keV. 
The boson number is N — 6 and S is an 0(6) label |17j . 



The relation of Eq. (|37[) ensures that sPo an d d m Po annihilate the states of the 0(6) ground 
band. The only three-body interactions that are partially solvable in 0(6) are thus PqUsPq and 

PqUiiPo- Since the combination Pj (n s + u^Pq = (N — 2)PqPo is completely solvable in 0(6), 
we can transcribe the 0(6)-PDS Hamiltonian in the form 

Hpds = Hds + V PqKsPo , (41) 



Table 3. Observed (EXP) and calculated B(E2) values (in e 2 b 2 ) for 196 Pt. For both the exact (DS) 
and partial (PDS) 0(6) dynamical symmetry calculations, the E2 operator is T(E2) = es[n( 2 ) +\U^ ] 
with es = 0.151 eb and x — 0.29. Only the state 0^ has a mixed 0(6) character [17] . 
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Here the dynamical symmetry Hamiltonian, Hds, is that of Eq. (|34p . since no new terms are 
added to it at the level of three-body interactions. The solvable states are members of the 
7-unstable deformed ground band 



The experimental spectrum and E2 rates of 196 Pt are shown in Fig. [2] and Table El The 0(6)- 
DS limit is seen to provide a good description for properties of states in the ground band (S = N). 
This observation was the basis of the claim |28j that the 0(6)-DS is manifested empirically in 
196 Pt. However, the resulting fit to energies of excited bands is quite poor. The Of, O3", and 
O4" levels of 196 Pt at excitation energies 0, 1403, 1823 keV, respectively, are identified as the 
bandhead states of the ground (v = 0), first- (v = 1) and second- (v = 2) excited vibrational 
bands [28]. Their empirical anharmonicity, defined by the ratio R = E(y = 2)/E(v = 1) — 2, 
is found to be R = —0.70. In the 0(6)-DS limit these bandhead states have r = L = and 
S = N, N — 2, N — 4, respectively. The anharmonicity R = —2/(N + 1), as calculated from 
Eq. (|35p . is fixed by N. For N = 6, which is the appropriate boson number for 196 Pt, the 0(6)- 
DS value is R = —0.29, which is in marked disagreement with the empirical value. A detailed 
study of double-phonon excitations within the IBM, has concluded that large anharmonicities 
can be incorporated only by the inclusion of at least cubic terms in the Hamiltonian |29j . In 
the IBM there are 17 possible three-body interactions [Tj- One is thus confronted with the 
need to select suitable higher-order terms that can break the DS in excited bands but preserve 
it in the ground band. On the basis of the preceding discussion this can be accomplished by 
the 0(6)-PDS Hamiltonian of Eq. ()4ip . The spectrum of -HpDS is shown in Fig. [2j The states 
belonging to the £ = N = 6 multiplet remain solvable with energies (|42b]), which obey the same 
DS expression, Eq. (f35|) . States with £ < 6 are generally admixed but agree better with the 
data than in the DS calculation. For example, the bandhead states of the first- (second-) excited 
bands have the 0(6) decomposition £ = 4: 76.5% (19.6%), S = 2: 16.1% (18.4%), and £ = 0: 
7.4% (62.0%). Thus, although the ground band is pure, the excited bands exhibit strong 0(6) 
breaking. The calculated 0(6)-PDS anharmonicity for these bands is R = —0.63, much closer to 
the empirical value, R = —0.70. It should be emphasized that not only the energies but also the 
wave functions of the £ = iV states remain unchanged when the Hamiltonian is generalized from 
DS to PDS. Consequently, the E2 rates for transitions among this class of states are the same 
in the DS and PDS calculations. Thus, the additional three-body term in the Hamiltonian (|4ip . 
does not spoil the good 0(6)-DS description for this segment of the spectrum. This is evident 
in Table [3] where most of the E2 data concern transitions between £ = iV = 6 states. 

5. Coexistence of partial dynamical symmetries 

The examples considered in previous sections involved Hamiltonians with a single PDS, 
describing stable structures, e.g., well-deformed nuclei. Multiple partial dynamical symmetries 
can occur in systems undergoing quantum phase transitions (QPTs) [30J. The latter are 
structural changes induced by a variation of parameters in the Hamiltonian. Such ground-state 
phase transitions are a pervasive phenomenon observed in many branches of physics [3T], and are 
realized empirically in nuclei as transitions between different shapes. In the IBM, the dynamical 
symmetry limits correspond to possible phases of the system and the relevant Hamiltonians for 
studying shape-phase transitions involve terms with from different DS chains [15j . The nature of 
the phase transition is governed by the topology of the surface En(/3, 7), Eq. ([TJ, which serves as 
a Landau's potential with the equilibrium deformations as order parameters. The surface at the 
critical-point of a first-order transition is required to have two degenerate minima, corresponding 
to the two coexisting phases. Specifically, the first-order critical surface is 



\[N](N)(T)n A L) , t = 0,1,2,..., N 
Epds = Bt(t + 3) + CL(L + 1) . 



(42a) 
(42b) 



E N (P, 1 = 0) 



2h 2 N(N - 1)(1 + /3 2 )" 2 /? 2 (/3 - A) ? 



(43) 
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Figure 3. Spectrum of H cr i(l3 = y/2), Eq. flU]), 
with h 2 = 0.05, N = 10. L(K = 0i) and L(K = 
2i) are the solvable SU(3) states of Eqs. (PH - flU 
with k — and fc = 1, respectively. L — O2, 3i are 
the solvable U(5) states of Eq. ® [30]. 
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Figure 4. U(5) (n d ) and SU(3) [(A,/*)] 
decomposition for selected spherical and de- 
formed states in Fig. [3] [30] . 



and has degenerate spherical and deformed minima at (3 = and (/3 = /?o>7 = 0), corresponding 
to spherical and axially-deformed shapes. A barrier of height h = h 2 N(N — 1)(1 — yl + /3 2 ) 2 /2 
separates the two minima. Such a surface can be obtained from Eq. (UJ with the following 
critical-point Hamiltonian 



H C ri((3o) 



Wo) , PLiPo) = A) + V7 (S S)% 



(44) 



P 2 m(Po) annihilates the condensate of Eq. { 
L, I /?o; iV, L), projected from it 

P 2m (/3o)|/3o,7 = 0;iV) 
P 2 m{/3o)\Po;N,L) 



with (/3 = /3q , 7 = 0) as well as the states of good 



, 
, 



L = 0,2,4,... ,2N . 



(45a) 
(45b) 



Consequently, H C ri(Po) nas a solvable zero-energy prolate-deformed ground band, composed of 
the L-projected states mentioned above 



\Po;N,L) E 
The multipole form of Hcri(fio) 



L = 0,2,4,... ,2N 

is given by 



(46) 



ho 



2(/3 2 A> - 2)h d + 2(1 - /3 2 )n 2 + 2(72(0(5)) - C 2 (0(3)) + Vup Il {2) • ^ {2) 1 , (47) 



where the various operators are defined in Table [TJ The fi d , n% C 2 (0(5)) and 6*2(0(3)) terms 
in Eq. (|47|) belong to the dynamical symmetry Hamiltonian of the U(5) chain 



U(6) D U(5) D 0(5) D 0(3) 

■i - "T" 4" "T 

[iV] (n d ) (r) n A L 



(48) 



and describe the dynamics of an anharmonic spherical vibrator. The II^ 2 ) • term can connect 
states with An d = ±1 and At = ±1,±3, hence breaks the U(5) DS. Nevertheless, H cr i((3o) has 
selected solvable eigenstates with good U(5) symmetry, 



\N,n d = t = L = 0) E = 0, 

\N,n d = r = L = 3) E = 6h 2 [${N-3) + 5\ , 



(49a) 
(49b) 
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tion for selected spherical and deformed states 
in Fig. [5] [30]. 



Figure 5. Spectrum of I? c h(A) = 1), Eq. ([H]). 
with ft 2 = 0.05 and N = 10. = 0i) are 

the solvable states of Eq. (|46|) with /3q = 1, which 
have good 0(6) symmetry a — N, but broken 0(5) 
symmetry. L = 02,3i are the solvable U(5) states 
of Eq. (01]) [3D]. 
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Figure 7. 0(5) (r) decomposition for the L = 0, 2 
states, Eq. (|4"6"]) with /3o = 1, members of the ground 
band (A" = 0i) of H cr i{Po = !)• Both states have 
0(6) symmetry a = N [3D]. 



and therefore, by construction, H cr i(/3o) exhibits U(5)-PDS of type I. 

For = y/2, the critical Hamiltonian of Eq. (|44p is recognized to be a special case of the 
Hamiltonian of Eq. (|18p . shown to have SU(3)-PDS of type I. As such, it has a subset of solvable 
eigenstates, Eqs. (|19p-(|20p. which are members of deformed ground g(K = 0) and "f k (K = 2k) 
bands with good SU(3) symmetry, (A,/x) = (2N — 4k, 2k). In addition, H cr i({3o = y/2) has the 
spherical states of Eq. (|49p. with good U(5) symmetry, as eigenstates. The spherical L = state, 
Eq. (|49a|) . is exactly degenerate with the SU(3) ground band, Eq. (fTU]l . and the spherical L = 3 
state, Eq. (|49bj) . is degenerate with the SU(3) 7-band, Eq. (|2Up with k = 1. The remaining 
levels of H cr i((3o = V2), shown in Fig. [3] are calculated numerically and their wave functions are 
spread over many U(5) and SU(3) irreps, as is evident from Fig. 2] This situation, where some 
states are solvable with good U(5) symmetry, some are solvable with good SU(3) symmetry and 
all other states are mixed with respect to both U(5) and SU(3), defines a U(5) PDS of type I 
coexisting with a SU(3) PDS of type I. 

For fa = 1, the critical Hamiltonian of Eq. (|47p involves the Casimir operators of 0(5) and 
0(3) which are diagonal in the corresponding quantum numbers, a and t, of the 0(6)-DS chain, 
Eq. (|33p . It also contains a term involving fid which is a scalar under 0(5) but can connect 
states differing by Act = 0, ±2 and a II^ 2 ) • term which induces both 0(6) and 0(5) mixing 
subject to Act = 0, ±2 and At = ±1,±3. Although H cr i(f3o = 1) is not invariant under 0(6), 
it has a solvable prolate-deformed ground band, Eq. (|4"6"]l with /3q = 1, which has good 0(6) 



symmetry, (a) = {N), but broken 0(5) symmetry. In addition, H cr i(f3o = 1) has the spherical 
states of Eq. (j4"9j) . with good U(5) symmetry, as eigenstates. The remaining eigenstates in Fig. [5] 
are mixed with respect to both U(5) and 0(6), as is evident from the decomposition shown in 
Fig. El Apart from the solvable U(5) states of Eq. (j4"9|) . all eigenstates of H cr i((3o = 1) are mixed 
with respect to 0(5) [including the solvable 0(6) states of Eq. (|46p with (3q = 1, as shown in 
Fig. [7] . It follows that the Hamiltonian has a subset of states with good U(5) symmetry and a 
subset of states with good 0(6) but broken 0(5) symmetry, and all other states are mixed with 
respect to both U(5) and 0(6). These are precisely the required features of U(5) PDS of type I 
coexisting with 0(6) PDS of type III. 

Second-order quantum phase transitions between spherical and deformed 7-unstable nuclei, 
can be accommodated in the IBM, by mixing the DS Hamiltonians of the U(5) chain, Eq. (|48|) . 
and the 0(6) chain, Eq. ([53"j) , Since U(5) and 0(6) are incompatible, yet both chains have a 
common 0(5) D 0(3) segment, we encounter a situation similar to that described in Eq. (0), 
which gives rise to 0(5)-PDS of type II [Tg] . 

6. Concluding remarks 

The notion of partial dynamical symmetry generalizes the concepts of exact and dynamical 
symmetries. In making the transition from an exact to a dynamical symmetry, states which are 
degenerate in the former scheme are split but not mixed in the latter, and the block structure 
of the Hamiltonian is retained. Proceeding further to partial symmetry, some blocks or selected 
states in a block remain pure, while other states mix and lose the symmetry character. A partial 
dynamical symmetry lifts the remaining degeneracies, but preserves the symmetry-purity of the 
selected states. 

Having at hand concrete algorithms for identifying and constructing Hamiltonians with PDS, 
is a valuable asset. It provides selection criteria for the a priori huge number of possible 
symmetry-breaking terms, accompanied by a rapid proliferation of free-parameters. This is 
particularly important in complicated environments when many degrees of freedom take part 
in the dynamics and upon inclusion of higher-order terms in the Hamiltonian. Futhermore, 
Hamiltonians with PDS break the dynamical symmetry (DS) but retain selected solvable 
eigenstates with good symmetry. The advantage of using interactions with a PDS is that they 
can be introduced, in a controlled manner, without destroying results previously obtained with 
a DS for a segment of the spectrum. These virtues greatly enhance the scope of applications of 
algebraic modeling of quantum many-body systems. 

PDSs appear to be a common feature in algebraic descriptions of dynamical systems. They 
are not restricted to a specific model but can be applied to any quantal systems of interacting 
particles, bosons, as demonstrated in the present contribution, and fermions [13 |, l32| . r33] . They 
are also relevant to the study of mixed systems with coexisting regularity and chaos, where they 
lead to a suppression of chaos [M1I35] . 
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